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Abstract. The Pliocene Model Intercomparison Project 
(PlioMIP) is the first coordinated climate model comparison 
for a warmer palaeoclimate with atmospheric CO 2 signifi- 
cantly higher than pre-industrial concentrations. The simula- 
tions of the mid-Pliocene warm period show global warming 
of between 1.8 and 3.6 °C above pre-industrial surface air 
temperatures, with significant polar amplification. Here we 
perform energy balance calculations on all eight of the cou- 
pled ocean-atmosphere simulations within PlioMIP Experi- 
ment 2 to evaluate the causes of the increased temperatures 
and differences between the models. In the tropics simulated 
warming is dominated by greenhouse gas increases, with the 
cloud component of planetary albedo enhancing the warm- 
ing in most of the models, but by widely varying amounts. 
The responses to mid-Pliocene climate forcing in the North- 
ern Hemisphere midlatitudes are substantially different be- 
tween the climate models, with the only consistent response 
being a warming due to increased greenhouse gases. In the 
high latitudes all the energy balance components become 


important, but the dominant warming influence comes from 
the clear sky albedo, only partially offset by the increases in 
the cooling impact of cloud albedo. This demonstrates the 
importance of specified ice sheet and high latitude vegeta- 
tion boundary conditions and simulated sea ice and snow 
albedo feedbacks. The largest components in the overall un- 
certainty are associated with clouds in the tropics and polar 
clear sky albedo, particularly in sea ice regions. These simu- 
lations show that albedo feedbacks, particularly those of sea 
ice and ice sheets, provide the most significant enhancements 
to high latitude warming in the Pliocene. 


1 Introduction 

Atmospheric carbon dioxide concentrations continue to rise 

due to anthropogenic emissions. The latest measurements 

show that annual mean concentrations have risen beyond 

390 parts per million (Conway et al., 2012). The Pliocene 
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Table 1 . Key model and experimental design parameters for each of the eight PlioMIP Experiment 2 simulations. 


GCM 

Atmospheric 

resolution 

(° lat x 0 long x levels) 

Ocean resolution 
(° lat x 0 long x levels) 

Boundary 

conditions 

employed 

Ocean 

initialization 

Reference 

CCSM4 

0.9 x 1.25 x 26 

1 x 1 x 60 

Alternate 

PRISM3 

(anomaly) 

Rosenbloom et 
al. (2013) 

COSMOS 

3.75 x 3.75 x 19 

3 x 1.8 x 40 

Preferred 

PRISM3 

(anomaly) 

Stepanek and 
Lohmann (2012) 

GISS-E2-R 

2 x 2.5 x 40 

1x1.25x32 

Preferred 

PRISM3 

Chandler et al. 
(2013) 

HadCM3 

2.5 x 3.75 x 19 

1.25 x 1.25 x 20 

Alternate 

PRISM2 
mPWP control 

Bragg et al. 
(2012) 

IPSLCM5A 

3.75 x 1.9x39 

0.5-2x2x31 

Alternate 

Pre-industrial 

control 

Contoux et al. 
(2012) 

MIROC4m 

2.8x2.8x20 

0.5-1. 4 x 1.4x43 

Preferred 

PRISM3 

Chan et al. 
(2011) 

MRI-CGCM 2.3 

2.8x2.8x30 

0.5-2 x 2.5 x 23 

Alternate 

PRISM3 

(anomaly) 

Kamae and 
Ueda (2012) 

NorESM-L 

3.75 x 3.75 x2 6 

3 x 3 x 30 

Alternate 

Levitus 

Zhang et al. 
(2012) 


was the last period of Earth history with similar to modern at- 
mospheric CO 2 concentrations (Kurschner et al., 1996; Seki 
et al., 2010; Pagani et al., 2010; Bartoli et al., 2011). These 
were associated with elevated global temperatures in both the 
ocean (Dowsett et al., 2012) and on land (Salzmann et al., 
2013). As the last period of global warmth before the climate 
transition into the bipolar ice age cycles of the Pleistocene, 
the mid-Pliocene warm period (mPWP) has been a target 
for both palaeoenvironmental data acquisition and palaeo- 
climate modelling over a number of years (Dowsett et al., 
1992; Chandler et al., 1994; Haywood et al., 2009; Dowsett 
et al., 2010). Although a number of different general circula- 
tion models (GCMs) have been used to simulate Pliocene cli- 
mates (Chandler et al., 1994; Sloan et al., 1996; Haywood et 
al., 2000, 2009; Dowsett et al., 201 1), it is only recently that a 
coordinated multi-model experiment has been initiated, with 
standardized design for mid-Pliocene simulations (Haywood 
et al., 2010, 2011). 

The Pliocene Model Intercomparison Project (PlioMIP) 
represents the first coordinated multi-model experiment to 
simulate a warmer than modem palaeoclimate, with high at- 
mospheric CO 2 concentrations (405ppmv). It has recently 
been added to the Paleoclimate Model Intercomparison 
Project (PMIP; Hill et al., 2012) and the first phase, in- 
corporating two simulations, completed. This paper focuses 
on PlioMIP Experiment 2, designed for coupled ocean- 
atmosphere GCMs (Haywood et al., 2011). Although, many 
of the large-scale features of the simulated Pliocene climate 
have been well documented (Dowsett et al., 2012; Haywood 
et al., 2013; Zhang et al., 2013, 2013a, b; Salzmann et al., 
2013), the causes of the simulated changes and differences 


between the simulations have not been extensively explored 
prior to this study. In this paper the energy balance of the 
PlioMIP Experiment 2 simulations are analysed in order to 
understand the causes of Pliocene atmospheric warming and 
the latitudinal distribution of increased surface air tempera- 
tures. This analysis allows us to analyse the causes of the 
warming, both directly through the simulated energy balance 
components and through examination of Earth system com- 
ponent changes that are driving these. It is also important 
when discrepancies with available proxy reconstmctions of 
Pliocene warming are considered (Salzmann et al., 2013). 

2 Participating models 

Eight different modelling groups have submitted simulations 
to PlioMIP Experiment 2. All of these models are coupled 
ocean-atmosphere GCMs, but range in complexity and spa- 
tial resolution. Table 1 contains the details of each of the 
models’ simulation, including the resolution at which it was 
run, the boundary conditions employed and the model initial- 
ization. Each of the simulations is documented in much more 
detail in a separate paper within a special issue of Geoscien- 
tific Model Development , referenced in Table 1 . The general 
climate sensitivity of the model and the annual mean global 
warming produced in its PlioMIP Experiment 2 simulation is 
detailed in Table 2. Further details about the models can also 
be found in Haywood et al. (2013) and the references therein. 
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Table 2. The climate sensitivity and global mean annual surface 
air temperature warming of each of the models with simulations 
in the PlioMIP Experiment 2 ensemble. Climate sensitivity is the 
equilibrium global warming for a doubling of atmospheric CO 2 - 
the values quoted here are a general value for each model and do 
not refer to the particular set-up and initialization procedures used 
for PlioMIP 


GCM 

Climate 

sensitivity 

(°C) 

Mean annual 
mPWP SAT 
warming (°C) 

CCSM4 

3.2 

1.86 

COSMOS 

4.1 

3.60 

GISS-E2-R 

2.7 

2.24 

HadCM3 

3.1 

3.27 

IPSFCM5A 

3.4 

2.03 

MIROC4m 

4.1 

3.46 

MRI-CGCM 2.3 

3.2 

1.84 

NorESM-F 

3.1 

3.27 


3 PlioMIP Experiment 2 

PlioMIP uses the latest iteration of the PRISM (Pliocene Re- 
search, Interpretation and Synoptic Mapping) mid-Pliocene 
palaeoenvironmental reconstruction, PRISM3 (Dowsett et 
al., 2010), as the basis for the imposed model boundary 
conditions. This reconstruction represents the peak averaged 
(Dowsett and Poore, 1991) warm climate of the mid-Pliocene 
warm period (mPWP; 3.246-3.025 Ma; Dowsett et al., 2010) 
in the middle of the Piacenzian Stage. It incorporates sea 
surface temperatures, bottom water temperatures (Dowsett 
et al., 2009), vegetation (Salzmann et al., 2008), ice sheets 
(Hill et al., 2007, 2010), orography (Sohl et al., 2009) and a 
global land-sea mask equivalent to 25 m of sea level rise. The 
vegetation, ice sheets and orographic reconstructions are all 
required as boundary conditions within the models, although 
they must be translated onto the resolution of each individ- 
ual model. Vegetation was reconstructed using the BIOME4 
classification scheme (Kaplan, 2001) and must therefore be 
translated onto the vegetation scheme used by each model. 

Although as part of PlioMIP a standard experimental de- 
sign was implemented, it was appreciated that not all of the 
modelling groups would be able to perform the ideal mPWP 
experiment. As such, alternate boundary conditions were 
specified for those models that could not effectively change 
the land-sea mask from the present-day configuration. This 
meant that the ocean advance specified in low-lying coastal 
regions and West Antarctica as well as the filling of Hudson 
Bay were not included in some of the simulations (Table 1). 
Furthermore a choice was given concerning the initial state 
of the ocean between a specification of the PRISM3 three- 
dimensional ocean temperatures (Dowsett et al., 2009) and 
initialization with the same ocean temperatures as the pre- 
industrial control simulation (Haywood et al., 2011). 


4 PlioMIP Experiment 2 global warming 

Overall the PlioMIP models simulate mPWP annual mean 
global surface air temperature (SAT) increases of 1.8-3. 6 °C 
(Table 2). Tropical temperatures increased by only 1.0- 
3.1 °C, while in the Arctic surface air temperatures increased 
by 3.5-13.2 °C (Fig. 1). Sea surface temperatures (SSTs) fol- 
low a similar pattern, but with a reduced magnitude of global 
warming and significantly greater warming in the North Pa- 
cific (Fig. Id). The patterns of warming in the northern mid- 
latitudes and southern high latitudes are much more variable 
between the different models. Relative variation between the 
models is largest in the North Atlantic, midlatitude mountain 
regions and central Antarctica for SATs (Fig. lc) and in the 
North Atlantic, North Pacific and sea ice areas of the Arctic 
and Southern oceans for SSTs (Fig. If). 

The warming of the PlioMIP simulations is accompanied 
by increased precipitation (Haywood et al., 2013) and mon- 
soonal activity (Zhang et al., 2013) and reductions in sea 
ice (Clark et al., 2013), although the Atlantic Meridional 
Overturning Circulation shows little response (Zhang et al., 
2013b). Global mean temperature response (Table 2), as well 
as polar amplification (Salzmann et al., 2013), do not show 
a strong correlation to either the use of preferred or alter- 
nate boundary conditions or to the initial conditions of the 
ocean. Although land-sea masks vary between the different 
models in the Hudson Bay and West Antarctic region, they 
do not show the largest relative variance, suggesting that the 
alternatives used in PlioMIP Experiment 2 do not introduce 
significant biases. 

5 Energy balance approach 

Energy balance analyses have been used in many palaeocli- 
mate simulations and ensembles to understand the simulated 
temperature changes (e.g. Donnadieu et al., 2006; Murakami 
et al., 2008). The results from each of the GCMs can be bro- 
ken down in to the various components in the energy balance 
of each individual simulation. The approach taken builds on 
the energy balance modelling of Heinemann et al. (2009) 
and Lunt et al. (2012), where globally averaged temperatures 
are approximated using planetary albedo a and the effective 
longwave emissivity s. 

^(1 - a) = eoT 4 , (1) 

where So is the total solar irradiance (1367Wm -2 ) and a 
is the Stefan-Boltzmann constant (5.67 x 10 -8 Wm -2 K -4 ). 
Planetary albedo is the ratio of outgoing (f) to incoming (|) 
shortwave radiation at the top of the atmosphere (TOA) and 
effective longwave emissivity the ratio of TOA to surface 
(SURF) upward longwave radiation, 

S\yj OA LW| 

a = s = TQA (2) 

SWjoa ' LW sV 
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Zonal mean Pliocene warming (°C) 
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Fig. 1. Multi-model mean PlioMIP Experiment 2 warming between mid-Pliocene and pre-industrial simulations, (a) Annual mean surface 
air temperature (SAT) warming, (b) zonal mean SAT warming (solid line), with shading showing the range of model simulations, and 
(c) relative variance between the PlioMIP Experiment 2 simulations (<r/ASAT). (d) Annual mean sea surface temperature (SST) warming, 
(e) zonal mean SST warming and (f) relative variance of SSTs. 


This can be expanded to approximate the one dimensional, 
zonally averaged temperatures at each latitude of the model 
grid by including a component for the implied net meridional 
heat transport divergence ( H ). 

SWjoa (1 -«) + // = so t\ (3) 

where 

H = - ((swj OA - SW| OA ) - LW| oa ) . (4) 


further broken down into the changes in the impact of at- 
mospheric greenhouse gases (Ar gge ), clouds (via impacts 
on both emissivity, A T cs , and albedo, A T ca ; see Sect. 6 
for discussion) and clear sky albedo (A T csa ; generally dom- 
inated by changes in surface albedo, but including atmo- 
spheric absorption and scattering components). In experi- 
ments and latitudes where changes in topography occur be- 
tween the Pliocene and pre-industrial times, the impact of 
these changes in surface altitude (ATfopo) must also be ac- 
counted for. 


Thus the temperature at each latitude in a GCM experiment 
is given by 


AT ^ A Tgg S + AT CS + A T ca + AT csa + ATjj ATtopo (T) 


T = 


SW 


TOA 


(1 - a) - H 


1/4 


EG 


= T(e, a, H). 


( 5 ) 


By applying the notation of Lunt et al. (2012) to denote the 
pre-industrial control experiment as a second experiment rep- 
resented by an apostrophe, the Pliocene surface air tempera- 
ture warming (AT) can be calculated by 

AT = T(s, a , H) - T (s' , a', H r ). (6) 


Due to their small changes relative to their absolute values, 
Pliocene warming can be approximated by a linear combi- 
nation of changes in emissivity (A T s ), albedo (A T a ) and 
heat transport (A7 h). However, these components can be 


Each of these components can be calculated from various 
combinations of Pliocene and pre-industrial albedos, emis- 
sivities and implied heat transports, although to differentiate 
between cloud and clear sky components some must be cal- 
culated in the clear sky case (denoted with a subscript cs). 


A7g g , = T(e cs , a, H) - T(4, a, H) - A T topo (8) 

A T cs = (T(s, a, H) — T (e cs , a, H)) 

-(T{s’,a,H)-T « s , «,//)) ( 9 ) 

A T ca = ( T(s , a, H) — T (s, a cs , H)) 

- (T(e, a', H) - T (e, a' cs , H)) (10) 

AT csa = T (e, « cs , H)-T (e, a' cs , H ) ( 1 1 ) 

A7h = T(s, a, H) - T(s, a, H') (12) 
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Changing the topography within a climate model can have 
many effects on the energy balance of the model, including 
changing circulation patterns, heat transport, surface condi- 
tions, cloud formation, storm generation, etc. Most of these 
features cannot be properly quantified without performing 
further simulations. However, a simple lapse rate correc- 
tion will remove the direct impact on surface temperatures 
of changing the height of the surface. Although lapse rates 
vary over time and space, the impact of changing the topog- 
raphy in the Pliocene simulations (Ar t0 po) can be approxi- 
mated by multiplying the change in topography (Ah) by a 
constant atmospheric lapse rate (y ~ 5.5 K km -1 ; Yang and 
Smith, 1985). 

ATtopo = Ah • y (13) 

6 Treatment of clouds within the energy balance 
calculations 

The energy balance calculations presented here split the plan- 
etary albedo and emissivity impacts into a component due 
to clouds and a clear sky one. In order to do this the clear 
sky radiation fluxes, which are the radiation fluxes that would 
have occurred had there not been any clouds, are used within 
the calculations. The cloud components are then the global 
temperature change due to the impact of clouds on planetary 
albedo and emissivity. If all else in the climate system re- 
mained the same except for the surface albedo then the clear 
sky albedo would show a change in temperature, but so also 
would the cloud albedo. This is because the impact of iden- 
tical clouds on the total planetary albedo have changed due 
to the surface which they are covering being different. Simi- 
larly an increase in highly reflective cloud types would have 
a greater impact on temperature over the dark oceans than 
over a high albedo land surface. 

A more thorough methodology for calculating the cloud 
and clear sky components would be to incorporate cloudi- 
ness into the calculations and thereby remove the reliance of 
cloud impact on the surface conditions. This requires an un- 
derstanding of the impact of different modelled cloud types 
on the energy balance and, as we are considering a multi- 
model ensemble, the consistent application of this across all 
the simulations within the ensemble. Total cloud fractions 
are produced by each model, but these use different overlap 
functions, all of which are simplifications of the observed de- 
pendencies of cloud overlap on atmospheric dynamics (Wang 
and Dessler, 2006; Naud et al., 2008). Therefore, the total 
cloud fractions of the models are not only incomparable, but 
may also not be the best representation of cloud impacts on 
energy balance. 

The impact of using the pragmatic methodology outlined 
above can be approximated using the surface downward 
shortwave radiation to approximate the impact of clouds. 
Figure 2 shows the albedo components of Pliocene warming 
in the HadCM3 model, using both this approximation and the 



Fig. 2. Comparison of the impact on Pliocene warming of the com- 
ponents of planetary albedo under different formulations for the 
HadCM3 model. The grey lines are the formulations presented in 
this Sect. 5, while the black lines represent an altered calculation 
of the albedo components approximating the impact of clouds us- 
ing the difference between the incoming radiation at the top of the 
atmosphere and the modelled downward shortwave radiation at the 
surface. 

methodology outlined in Sect. 5. Although using the outlined 
energy balance approach produces greater magnitudes of im- 
pacts from clear sky and cloud albedo, as would be expected, 
the overall structure and the conclusions that would be drawn 
from these calculations remain. 


7 Energy balance results for individual simulations 

The energy balance calculations for each of the individual 
simulations within the PlioMIP Experiment 2 ensemble are 
shown in Fig. 3. The overall structure of the energy balance 
components is largely the same between all the simulations, 
although there are differences in the responses in each of 
the models and their magnitudes. In the tropics the warm- 
ing is dominated by the greenhouse gas emissivity, with the 
other components having a small impact on warming. Those 
simulations showing greater tropical warming tend to also 
have a significant warming component from cloud impact 
on albedo. The midlatitudes are the region where the mod- 
els show the least consistency, especially in the Northern 
Hemisphere. The only consistent component of the warming 
comes from greenhouse gas emissivity, although the cloud 
and clear sky albedos also tend to have a warming impact. In 
the high latitudes much of the warming comes from changes 
in the clear sky albedo, which are only partially offset by 
changes in the overall cloudy sky impact on planetary albedo. 
There are also slightly enhanced greenhouse gas warming 
and a warming impact of cloud emissivity. In the Southern 
Hemisphere the clear sky albedo warming often has a dou- 
ble peak, the first around 60° S, representing changes in the 
simulated Southern Ocean sea ice and the second over the 
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Fig. 3. Energy balance analysis for each of the eight PlioMIP Experiment 2 simulations, from (a) CCSM4, (b) COSMOS, (c) GISS-E2-R, 
(d) HadCM3, (e) IPSLCM5A, (f) MIROC4m, (g) MRI-CGCM 2.3 and (h) NorESM-L. Plots show the zonal mean warming, at each latitude 
in the model, from each of the energy balance components. The solid black line is the zonal mean surface air temperature increase from the 
GCM simulation, while the dashed grey line is the Pliocene warming approximated by the energy balance calculations. 
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Antarctic continents, from prescribed changes to the Antarc- 
tic ice sheets. 


8 PlioMIP Experiment 2 energy balance 

In order to evaluate the simulation of warm climates of 
the Pliocene in general, a simple mean of the energy bal- 
ance components from each of the individual simulations 
within the PlioMIP Experiment ensemble has been per- 
formed. When combined with the range of values within the 
ensemble this allows an assessment of the general cause of 
warming within the PlioMIP simulations and the robustness 
of any conclusion that can be drawn. Figure 4 shows the 
ensemble mean of the various energy balance components 
along with the range from the eight simulations, while Fig. 5 
shows the individual energy balance components for each of 
the PlioMIP simulations. 

Clear sky albedo includes contributions from surface 
albedo changes and atmospheric absorption and scattering. 
The latter could become important, even in models with 
no mechanisms for changing atmospheric transparency, as 
atmospheric thickness can increase due to changes in sur- 
face altitude. In the PlioMIP simulations clear sky albedo 
shows little contribution to warming in the tropics and South- 
ern Hemisphere midlatitudes. In the Northern Hemisphere 
midlatitudes most models show a warming due to clear 
sky albedo, apart from the MRI-CGCM 2.3 simulation that 
shows a cooling (Fig. 5a). In the polar regions, all the simu- 
lations show a strong warming signal from clear sky albedo, 
although the range in the magnitude of this warming is 
large. Changes in clear sky albedo mostly reflect changes on 
Earth’s surface. Vegetation, snow and ice (both terrestrial ice 
masses and sea ice) are generally the main contributors to 
these changes. The warming found in the Northern Hemi- 
sphere, from 15-60° N is largely being driven by changes in 
the vegetation boundary conditions, particularly over the Sa- 
hara, Arabia and central Asia (Fig. 6). In the Arctic, warm- 
ing due to clear sky albedo is primary driven by changes in 
ice sheet boundary conditions (reduced Greenland Ice Sheet) 
and changes in the predicted sea ice, but also by the pole- 
ward shift of the Arctic tree line (Salzmann et al., 2008). 
In the Southern Ocean and Antarctica the warming due to 
clear sky albedo has a double peak in most models, reflect- 
ing a reduction in the simulated Southern Ocean sea ice and 
a reduction in the prescribed Antarctic Ice Sheet (Fig. 7). 
Most models show reductions in sea ice coverage across the 
whole range of latitudes in which sea ice is modelled, except 
for the highest Northern Hemisphere latitudes and in CCSM 
close to Antarctica (Fig. 7c). The simulated changes in total 
cloud are relatively small and inconsistent between the mod- 
els (Fig. 7e), so do not appear to introduce systematic biases. 

All the simulations show a warming due to greenhouse gas 
emissivity of around 1-2 °C. These impacts are largely con- 
stant across latitudes, but with a slight polar amplification, 



Fig. 4. Summary of the PlioMIP Experiment 2 energy balance anal- 
ysis. Solid line shows the multi-model mean warming for each com- 
ponent, with the associated shading representing the range. These 
values have been interpolated onto a 1° latitude grid for compari- 
son purposes. 


especially in the Arctic (Fig. 5b). This is consistent with 
the prescribed increases in CO 2 (at 405 ppm for the mid- 
Pliocene, as opposed to 280 ppm in the pre-industrial sim- 
ulations). The amplified high-latitude response is due to in- 
creases in the atmospheric water vapour predicted by the 
models. Differences in the simulation of this water vapour 
increase between different models explain why the range 
of temperature increases due to greenhouse gas warming is 
much higher in the polar regions. 

The changes in the impact of clouds on planetary albedo 
are small in the tropics and midlatitudes. Different models 
seem to produce significantly different responses making the 
signal particularly noisy (Fig. 5c). However, the multi-model 
mean cloud albedo impact on warming appears to reflect 
some of the large-scale features of the PlioMIP simulations 
(Haywood et al., 2013). Between the Equator and ~ 45° there 
is a general warming due to a reduction in the impact of cloud 
albedo, interrupted by a cooling in the Northern Hemisphere 
tropics. This cooling is due to an increase in cloud cover re- 
sulting from a northward shift of the Inter-Tropical Conver- 
gence Zone (Kamae et al., 2011). In the high latitudes the 
increased impact of clouds on planetary albedo, partly due to 
changes in the underlying surface albedo, leads to a signif- 
icant cooling, peaking at between 3 and 6 °C in both hemi- 
spheres. Cloud emissivity shows a similar pattern of impacts, 
but in the opposite direction. However, the response is gen- 
erally of a smaller magnitude (Fig. 5d). 

Reconstruction of mid-Pliocene sea surface temperatures 
has led to increased heat transport in the North Atlantic being 
suggested as a primary driver of warming in the mid-Pliocene 
(Dowsett et al., 1992; Raymo et al., 1996). However, the im- 
plied overall meridional heat transport in the PlioMIP simu- 
lations, which integrates both oceanic and atmospheric trans- 
ports, shows little coherent signal. The fact that there is only 
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Fig. 5. Breakdown of the energy balance components, (a) clear sky albedo, (b) greenhouse gas emissivity, (c) cloud albedo, (d) cloud 
emissivity and (e) implied meridional heat transport. Solid black line shows the multi-model mean, range is shown by the grey shading and 
each of the individual model results is shown by the coloured solid lines. 


one region where all of the simulations show a temperature 
change of the same direction suggests that the only robust 
conclusion that can be drawn about heat transport is a re- 
duction of overall transport into the Arctic (Fig. 5e). This 
would be an expected result of reduced thermal gradients 
due to polar amplification in the Arctic region under climate 
warming. These energy balance calculations support analy- 
sis of the Atlantic Meridional Overturning Circulation in the 
PlioMIP ensemble, which shows that there is little change in 
the northward heat transport in the North Atlantic (Zhang et 
al., 2013b). This calls into question the role of ocean heat 
transport in the general warming of the mid-Pliocene. How- 
ever, it may be important in the Pliocene variability of sea 
surface temperatures, which is particularly high in the North 
Atlantic (Dowsett et al., 2012). 


9 Conclusions 

The mid-Pliocene was probably the last time in Earth’s his- 
tory when atmospheric carbon dioxide concentrations were 
similar to today (Kiirschner et al., 1996; Seki et al., 2010; 
Pagani et al., 2010; Bartoli et al., 201 1). It has been the focus 
of palaeoenvironmental reconstructions and palaeoclimate 
model experiments for many years. However, the recently 
begun Pliocene Model Intercomparison Project is the first 
time that coordinated multi-model experiments, with com- 
mon boundary conditions and experimental protocols, have 
been undertaken. The warming seen in the Pliocene has been 
well documented from a wide variety of sites from across 
the globe and using a number of different proxy techniques 
(Dowsett et al., 2012; Salzmann et al., 2013). Previous sim- 
ulations of Pliocene warmth have been performed with only 
a single model and multi-model analyses have been severely 
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Fig. 6. Spatial distribution of multi-model mean changes in the clear sky albedo between mid-Pliocene and pre-industrial simulations. This 
primarily shows changes due to specified vegetation and global ice sheets and modelled sea ice and snow cover, although it includes an 
atmospheric component. Greyscale shows reductions in albedo, generally associated with warming and yellow shading indicates increases 
in albedo that would generally cause a reduction in surface air temperature. Increases in albedo closely follow the northward expansion of 
grassland imposed on the simulations (Haywood et al., 2010). 


hampered by differing experimental designs (Haywood et al., 
2009). For the first time a robust analysis of the causes of 
warming in Pliocene climate models is possible. 

Energy balance calculations show that the tropical warm- 
ing seen in all the models is primarily caused by green- 
house gas emissivity, with specified increases in atmospheric 
CO 2 concentration being the most important factor. Along 
with different sensitivity to the imposed CO 2 concentrations, 
changes in warming due to the cloud impact on planetary 
albedo drive differences between the models in the tropics. 
At polar latitudes all the energy balance components become 
important, but clear sky albedo is the dominant driver of 
the high levels of warming and polar amplification. This is 
largely due to reductions in the specified ice sheets and sim- 
ulated sea ice, but in the Northern Hemisphere also reflects 
a northward shift in the treeline. The models show a very 
different response in the midlatitudes of the Northern Hemi- 
sphere, with large uncertainties in the relative contributions 
of the different energy balance components. This is partic- 
ularly true for the North Atlantic and Kuroshio Current re- 
gions, where intermodel variability is highest and warming 
is simulated very differently (Haywood et al., 2013). A more 
complete picture of these currents, their strength and variabil- 
ity within the Pliocene, would enable a much better analysis 
of the skill of the models in these key regions. 

Atmospheric CO 2 concentrations remain a significant un- 
certainty in Pliocene climate, with different proxy techniques 
producing values between pre- industrial levels and double 
pre-industrial levels (Kiirschner et al., 1996; Seki et al., 2010; 
Pagani et al., 2010; Bartoli et al., 2011). As tropical warming 
is largely driven by this factor, simulations with accurate rep- 
resentation of low latitude clouds could provide some new 


insight into the levels of CO 2 required to produce Pliocene 
climates (for example in a similar way to Lunt et al., 2012, for 
the Eocene). Alternatively, accurate reconstructions of sur- 
face temperatures and atmospheric CO 2 in combination with 
modelling studies could reveal the extent of changes to trop- 
ical cloud cover in the warmer Pliocene world. 

Particularly strong warming in the high latitudes is driven 
by changes in albedo, especially from sea ice, ice sheets, 
snow cover and vegetation, which are only partially offset by 
extra cooling components from increased cloud albedo im- 
pacts. This is the region with the largest warming signal and 
also the largest uncertainties between the simulations. There- 
fore, improvements in the reconstruction of global ice cover 
and Arctic vegetation, along with improved data to evaluate 
the simulation of sea ice and high Arctic atmospheric and 
ocean temperatures, could significantly improve the simu- 
lations and allow much better constraints on total Pliocene 
warming. From the PlioMIP Experiment 2 simulations it ap- 
pears that higher CO 2 concentrations warmed the planet dur- 
ing the Pliocene and drove large surface albedo feedbacks in 
the high latitudes through changes in sea ice, vegetation and 
ice sheets. The latter two of these factors are important com- 
ponents of long-term Earth system sensitivity, further sup- 
porting a long-term response of CO 2 greater than conven- 
tional climate sensitivity (Lunt et al., 2010; Haywood et al., 
2013). 
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Fig. 7. Zonal mean of a series of factors that contribute to the clear sky albedo and its impact on Pliocene warming, (a) The total Pliocene 
warming and the clear sky albedo component of this (black and green respectively; see Fig. 4). (b) Total zonal surface albedo change in each 
of the 8 models, (c) Zonal mean sea ice percent coverage in the models, solid line is the Pliocene total and the dashed line is the change 
between Pliocene and pre-industrial experiments, (d) Imposed snow free albedo change, as calculated from the observed and reconstructed 
Pliocene mega-biomes, which include vegetation and ice sheet changes (Salzmann et al., 2008; Haywood et al., 2010). (e) Change in modelled 
total percent cloud cover for each model. These are calculated by each of the models using their own algorithm and cloud parameterizations, 
so are indicative rather than directly comparable. 
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